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LIST  OP  SYMBOLS 


A 


k 


K 

n 

P 

Ap 

P 

r 


a  field  point  where  the  induced  flow  field  of  the  pro¬ 
peller  is  to  be  evaluated  (see  Figure  l) 
a  source  point  on  the  propeller  disk  (see  Figure  l) 


complete  elliptic  integral  as  defined  by  equation  [173 

constant  defined  by  equation  [7] 

’  T 


thrust  coefficient  = 


i? 


V^S 


argument  of  complete  elliptic  integral  defined  by 
equation  [15] 

argument  of  complete  elliptic  Integral  defined  by 
equation  [19] 

complete  elliptic  integral  of  the  first  kind 
power  index  of  the  assumed  loading  variation,  equation  [5] 
perturbed  pressure,  i.e.,  pressure  at  any  point  minus  the 
pressure  in  the  uniform  stream 
pressure  jump  across  the  propeller  disk 
Legendre  function 
radial  coordinate 


=  +  1,2  (see  Figure  l) 

propeller  radius 
propeller  disk  area 

perturbation  velocity  in  the  radial  direction 


v  ,v  ,v  perturbation  velocity  in  the  x-,  y-,  and  z-  direction 
x  y  z 

V  uniform  stream  velocity 

x,y, z  cartesian  coordinates 

6  angular  coordinate  of  the  field  point 

angular  coordinate  of  the  source  point 
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6  cos  ( — )  (see  Figure  l) 

o 

p  density  of  fluid 

$  velocity  potential  exterior  to  the  propeller  slipstream 

a)  distance  between  the  source  point  A1  and  field  point  A 

(see  Figure  l) 

Prime  denotes  differentiation  with  respect  to  the  argument. 
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ABS TRACT 


The  flow  field  of  an  infinitely  bladed 
propeller  of  disk  loading  proportional  to  (r/R)n*  (l-r/R) 
in  an  incompressible,  inviscid  fluid  is  expressed  in 
terms  of  a  series  of  Legendre  functions.  Numerical 
calculations  have  been  carried  out  for  the  case  n  =  5. 
The  radial  velocity  at  r  =  R,  and  the  shape  of  the 
stream  surface  passing  through  the  periphery  of  the 
propeller  disk  are  determined  and  compared  to  that  of 
uniform  disk  loading.  ^ 


1'  4  kv 
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THE  PLOW  FIELD  OF  AN  INFINITELY  BLADED  PROPELLER 
WITH  RADIALLY  NON-UNIFORM  LOADING 


Ic  Introduction, 

The  marine  propulsion  system  consisting  of  a  propeller  en¬ 
closed  by  a  nozzle  has  found  many  useful  applications,  especially 
when  high  disk  loadings  are  required 0  Since  the  advent  of  the 
supercavitating  propellers,  it  has  been  suggested  by  Tulin  that 
a  nozzle-supercavitating  propeller  system  can  extend  to  lower 
speeds  the  range  of  speeds  in  which  a  supercavitating  propeller 
can  be  practically  operated.  In  order  to  make  a  rational  design 
of  such  a  system,  and  to  predict  performance  both  at  design  and 
other  forward  speeds,  a  knowledge  of  the  flow  field,  especially 
of  the  radial  component  of  velocity  near  the  tip  of  the  propeller, 
is  needed. 

The  flow  field  of  an  infinitely  bladed  propeller  in  an  in¬ 
compressible,  inviscid  fluid  may  be  calculated  using  either  one 
of  the  following  two  equivalent  methods:  (l)  by  solving  the 
exact  continuity  equation  together  with  boundary  condition  at  the 
propeller  disk  derived  from  the  linearized  Bernoulli's  equation, 
or  (2)  by  solving  the  Laplace  equation  for  the  perturbed  pressure 
p  derived  from  the  linearized  equation  of  motion,  together  with 
the  specified  pressure  jump  across  the  propeller  disk.  With  both 
of  these  methods,  the  propeller  is  considered  to  be  lightly  loaded 
otherwise  the  equations  concerned  cannot  be  linearized.  The  prob¬ 
lem  then  reduces  to  the  solution  of  the  Laplace  equation  with 
specified  boundary  condi tionsj  solutions  may  be  obtained  either 
by  separation  of  variables  or  through  a  suitable  distribution  of 
singularities.  For  propellers  with  uniform  loading,  the  solution 
has  been  given  by  Ktichemann  and  Weber  (1953)#  Meyerhoff  and  Finkel 
stein  (1958)  and  others.  Tachmindji  (1959)  has  given  the  axial 
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velocity  field  of  a  propeller  with  optimum  (non-uniform)  loading. 

In  the  solution  of  the  uniformly  loaded  disk  (cf.  Kiichemann 
and  Weber,  1953  )>  the  radial  component  of  the  induced  velocity  ex¬ 
hibits  a  logarithmic  singularity  at  the  periphery  of  the  propel¬ 
ler  disk.  Since  logarithmic  singularities  are  integrable,  the 
streamline  near  the  propeller  tip  is  continuous  but  with  infinite 
slope.  The  region  influenced  by  this  singular  radial  velocity  is 
very  small  compared  to  the  radius  (less  than  0.0001  R),  but  never¬ 
theless  troublesome.  In  the  belief  that  this  singular  behavior 
arises  out  of  the  discontinuous  loading  at  the  periphery,  it  was 
decided  to  investigate  loadings  which  decrease  to  zero  continuously 
at  the  periphery  of  the  propeller  disk. 

In  this  report,  starting  with  the  Laplace  equation  govern¬ 
ing  the  perturbation  pressure  p,  we  obtain  the  solution  to  the 
problem  In  which  the  propeller  loading  varies  as  (r/R)n(l-r/R), 
where  R  is  the  radius  of  the  propeller.  In  terms  of  a  series  of 
Legendre  functions.  Numerical  calculations  have  been  carried  out 
for  the  case  n  =  Values  for  the  radial  component  of  the  per¬ 
turbation  velocity  at  r  =  R  are  determined  and  are  compared  with 
those  obtained  for  the  case  of  uniform  loading.  The  shape  of  the 
streamline  passing  through  the  propeller  tip  has  also  been  deter¬ 
mined  and  compared  with  the  case  of  uniform  loading.  Although  the 
radial  components  of  the  perturbation  velocities  very  close  to  the 
propeller  tip  for  these  two  loadings  are  quite  different  in  charac¬ 
ter,  the  resulting  streamline  shapes  are  on  the  whole  quite  similar. 

II.  Velocity  Field. 

The  linearized  equations  of  motion  for  the  perturbed 
physical  variables,  pressure  and  velocity,  in  rectangular  coordinates 
as  shown  in  Figure  1  ares 
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dv  dv  dv 

where  X  is  the  force  per  unit  volume  which  is  imparted  to  the 

fluid  by  the  propeller  disk,  V  is  the  velocity  of  the  uniform 

stream,  and  v  ,  v  ,  and  v  are  perturbed  velocity  components, 
x  y  z 

Prom  these  equations,  one  easily  obtains  the  Laplace  equation  for 
the  pressure  (except  at  the  propeller  disk  where  a  jump  in  pres¬ 
sure  occurs): 

v2p  =  o.  [3] 

In  axisymmetric  flow,  as  in  the  present  case,  equation  [3]  may 
be  written  in  the  cylindrical  coordinate  system  as 


where  x  and  r  are  the  axial  and  radial  coordinates  respectively. 
The  specified  pressure  jump  at  the  propeller  disk,  i.e„,  for  x  *  0 
and  r  <  R,  is  taken  for  the  loadings  under  consideration  as* 


in  which  CT  is  the  thrust  coefficient  of  the  propeller.  The  solu¬ 
tion  to  equation  [4]  satisfying  the  boundary  condition  [5]  can  be 
obtained  by  a  distribution  of  doublets,  with  strength  equal  to 
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the  pressure  jump,  at  the  propeller  disk.  We  shall  use  subscript  1 
to  denote  the  source  point  A^  (x1  =  0,  r^,  0^)  on  the  propeller  disk. 
Because  of  axi symmetry,  we  consider  the  field  point  A  (x,r,0)  to  be 
in  the  xy  plane,  as  shown  in  Figure  1.  The  distance  between  A  and  A^ 
is  indicated  by  cd.  For  the  sake  of  convenience,  we  consider  all  dis¬ 
tances  are  non-dimensionalized  with  respect  to  the  radius  of  the 
propeller  R  in  the  mathematical  treatment.  But  the  numerical  re¬ 
sults,  as  presented  in  plots,  are  presented  in  the  actual  physical 
coordinate  system. 

The  pressure  p  is  given  by 

a 


p(x'r)  =  h\  rl  (l-ri}  SET  fa  dSj 

s  1 


[6] 


where 


c  =  ipv2cT  lq+2-)-(n+3-).  , 


and  S  denotes  the  surface  area  of  the  propeller  disk.  Equation 
[6]  may  be  written  as 


2  w 


p(x*r)  -  If  s^; 


n+1 


(l-rx  )drx 


'0 


be. 


CD 


[7] 


[8] 


The  0^  integral  in  [8]  may  be  written  in  terms  of  Legendre  functions, 
as  shown  by  Kellogg  (1953 ), 

,27 r 


d0. 


CD 


27 r 
r 


JQ 


P2(cos0o)  +  (~)  P4(cos0o) 

o  o 


-  (p-)6  p6(oo»«„)  + ... 

0 


,  for  ro  >  r1  ,  [9] 
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2  2  2  x 

where  r  =  r  +  x  ,  and  cos0  =  —  ,  see  Figure  1.  Hence  at  any 
o  or 

o 

point  such  that  rQ  >  1,  we  have 


/  N  C  d  1 

p(x’r>  =  2  ST  - 
1  0 


n+2 )  [  nH 


n+4  j in 


17^51  ^  P2(ooseo) 


if 

+  2-4(n+6Hn+7)  ^r"^  VCOS0o)  +  '*• 


From  the  first  equation  of  [1],  we  obtain  for  points  outside  the 
slipstream 

v  =  £_ 

x  pV  ° 


[10] 


[11] 


Outside  the  slipstream,  there  exists  a  velocity  potential  $  which 
may  be  obtained  as 


d>  =  _ 


v  dx. 
x 


[12] 


Noting  the  relations 

f  dr  dr 

!  o  _ o 

'  S3T 


[13] 


dcos© 


dcosQ 


we  obtain  from  equation  [12] 


A 

■/ 


dx 

pV 


-00 
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$  =  - 


^pVr 


1 

'n+2 )  (nH 


(i-)  P0(cose  ) 


+  2* 4-  (n+6^(n+7 j  (rj  "  ?k<‘CoaBo'>  +  *’* 


for  x  >  0,  r  >  1. 


[14] 


This  result  is  exactly  the  same  as  that  of  Korvin-Kroukovsky(l956) 
for  the  special  case  of  uniform  loading.  Now, 

/  \ 

vr(x,r)=  -  ^ 


- (n+2)(n+3)<  -j-  (n+2)(n+3)  ' 

o  L 

4 

+  2-4- (n+6)(n+7)  ^  P^(COS0O)  +  ••• 


2(n+4](n+5)^r J  p2^cos0o^ 


1-3 

2*4* (n+6) (n+7 


(?")  P2(COS0O) 


(zr)  p1(co86„)  +  .. 


for  rQ  >  1,  [15] 


where  the  prime  denotes  differentiation  with  respect  to  the  argument. 

Equation  [11],  together  with  [10],  give  the  axial  component,  and 

equation  [15]  gives  the  radial  component  of  the  perturbed  velocity 

outside  a  sphere  of  radius  1  due  to  a  propeller  with  loading  given 

by  equation  [5].  The  axial  component,  (v  -v  ),  is  antisym- 

x  =  o 

metric  and  the  radial  component,  v  ,  is  symmetric  with  respect  to 
the  plane  x  =  0. 
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We  shall  now  let  n  =  5  in  the  loading  expression  [5].  The 
resulting  pressure  jump  along  a  radius  of  the  propeller  disk  is 
given  in  Figure  2.  This  distribution  is  chosen  for  the  following 
reasons:  (l)  the  occurrence  of  the  maximum  loading  near  the  tip, 

r/R  =  0.833*  may  closely  approximate  the  loadings  on  certain  pro¬ 
pellers  enclosed  by  a  nozzle,  and  (2)  the  near  zero  loading  near 
the  center  may  approximate  the  presence  of  a  hub.  For  this  loading 
then,  the  radial  component  of  the  perturbed  velocity  at  r  -  1  is 

-  1 


P2m(°) 

TwJWzJ 


(-1 J"*1  (awi  )P2m(coseo)-  (-1 ) 


m  _1_ 
2m 


p'  (cose  ) 
2m'  o' 


[16] 


For  x  >  0.4,  this  series  converges  quite  rapidly,  terms  up  to  and 
including  Pio(cose0)  are  used  to  evaluate  vr(x, l).  At  x  =  0,  the 

Legendre  functions  are  expressible  in  terms  of  gamma  functions  of 
integer  argument,  which,  for  large  values  of  m,  may  be  approximated 
by  Stirling’s  formula.  Using  the  integral  method  of  summing  up 
the  monotonically  decreasing  terms  for  large  m,  the  value  of  vr(0,l) 


is  obtained  with  a  maximum  error  of  3  percent.  For  0  <  x  <  0.4, 
because  convergence  of  the  series  is  quite  slow,  terms  up  to  and 
including  Pig.fcosQ  ),  as  tabulated  by  Tallqvist  (1937)*  are  used 
in  the  summation.  Error  estimates  are  very  difficult  to  make.  In 
order  to  check  the  results  thus  obtained,  we  resort  to  a  direct 
numerical  integration  of  equation  [8], 

First,  according  to  Kellogg  (1953)*  the  G integral  in  [8] 
is  evaluated  as 


j 


4 


Nondimensional  Pressure  Jump  across  Propeller  Disk; 
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The  integral  in  equation  [19]  is  evaluated  numerically  using  Simp¬ 
son's  rule  for  points  0  <  x  <  0.4.  The  results  of  vr(x,  1  J/C^V  are 
within  10  percent  of  those  obtained  by  summing  the  series.  The 
values  of  v^(x, l)  are  plotted  in  Figure  3. 


The  flow  field  of  a  uniformly  loaded  propeller  has  been  given 
by  Kdchemann  and  Weber  (1953).  In  particular,  downstream  of  the 
propeller  v^(x,  l)  is  expressed  as 


Figure  3-  Axial  Variation  of  the  Induced  Radial  Velocity  at  Propeller  Radius 
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vr(x,l) 

CTV 


2 

(1  -  ~)  K(kx)  -  EfkpJ 


[21 J 


in  which  E  is  the  complete  elVptic  integral  of  the  second  kind. 


and 


k2  =  _t_ 

1  x2+4 


[22] 


It  Is  seen  that  as  x -*  0,  vr(x, l)  behaves  like  In  xj  it  becomes 

singular  at  x  -  0,  as  it  is  shown  in  Figure  3*  This  singular  be¬ 
havior  may  be  expected  for  loadings  which  decrease  to  zero  discon- 
tinuously  at  the  propeller  tip.  This  logarithmic  singularity,  as 
it  is  well  known,  is  integrable.  Therefore,  the  streamline  is  con¬ 
tinuous  at  the  propeller  tip,  but  with  an  infinite  slope.  In 
reality,  this  singular  behavior  will  not  be  present  since  the  load¬ 
ing  may  decrease  to  zero  quite  sharply,  but  never  discontinuously. 
Furthermore,  the  viscous  nature  of  the  fluid  would  prohibit  the 
presence  of  infinite  velocities.  It  is  seen  from  Figure  3>  that 
for  x/R  >  0.3,  the  values  of  vr(x, l)  due  to  these  two  different 
loadings  are  quite  similar. 

III.  Streamline  Shape. 


Along  a  streamline,  the  following  relation  holds: 

dr  _  vr 
dx  “  V+v 


[23] 


Consistent  with  the  assumptions  made  to  arrive  at  the  linearized 
equation  of  motion,  v  may  be  neglected  comparing  with  V.  The 
streamline  shape  is  then. 


r-1  =  /  7f-  dx. 


[24] 


* 
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Using  the  results  obtained  in  the  last  section  for  the  radially 
non-uniform  loading,  this  above  integration  is  performed  numer¬ 
ically  and  the  results  are  shown  in  Figure  4.  In  the  uniformly 

,  j  re¬ 
loaded  case,  we  expand  both  K  and  E  in  terms  of  =  V  1-k^  ; 


then  by  termwise  integration,  we  obtain 


r  -  1 


i  .3 


0.386  -  In  kj  k|  +  (0.04-0.75  In  k1)(k1)  +.  . . 


>[25] 


The  results  are  also  shown  in  Figure  4.  It  is  seen  that  the 
maximum  difference  between  these  two  results  is  about  0.01  C^. 

However,  this  may  incur  as  much  as  14  percent  relative  difference 
between  the  mean  radial  velocities  calculated  for  these  two 
different  loadings.  Both  loadings  lead,  of  course,  to  the  same 
asymptotic  value  of  r  at  distances  far  from  the  propeller. 
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